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Abstract: We study the localization of gravity on a string-like topological defect within a 
6-dimensional space-time. Assuming zero cosmological constant we find complete numeri- 
cal solutions to a set of first-order, Bogomol'nyi-Prasad-Sommeferld (BPS)-like, equations 
for the metric and the scalar field, where the dynamics of the latter are dictated by a 
super gravity- type potential. Our axially symmetric solutions have no deficit angle and 
factorize as AdS^ x Si far from the core. They are regular everywhere, providing complete 
smooth cigar-like geometries. The total energy of these configurations depends only on the 
boundary conditions for the warp factor and it is shown to vanish. 
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1. Introduction 

Already since the early eighties jl], |2|, ^ physicists have toyed with the idea of considering 
our 4D Universe as the internal space of a topological defect living in higher dimensional 
space-time. Such defect could either be a domain wall (if the total of space-time dimensions 
is equal to five), a global string (in six dimensions), a monopole (in seven dimensions) or an 
instanton (in eight dimensions). This idea received further boosts, first with the realization 
that solitons in string theory (i.e. D-branes) were ideal candidates to localize gauge and 
matter fields followed by the work of Randall and Sundrum Q on the localization of 
gravity in a three-brane domain wall in five space-time dimensions. 

Since then exhaustive work has been devoted to the study of five-dimensional models 
and slightly less to higher dimensional setups. In particular, for D=6 or the global string 
case, a number of issues have already been discussed. Whereas in D=4 the global string in 
the presence of a static metric is singular [^], and finite when we include time-dependence 
0, we see that the corresponding D=6 model continues to be singular when the metric is 
time-independent and in the presence of a zero cosmological constant, while it becomes 
regular again if a negative cosmological constant is introduced Q. So far we are always 
assuming a system with gravity and a scalar field with a U(l) invariant interaction. Let 
us mention that, already, a fair amount of work has been devoted to the study of the local 



string in D=6 [10| and, in particular, to finding numerical solutions to the Abelian Higgs 



model 



Also, a very important aspect of these D=6 models is their potential phenomenological 



applicability. In that respect the localization of matter fields has been studied in refs. |12] 



P]. Other issues, such as their cosmology or the status of the cosmological constant 
problem [17|-|23] have also been considered elsewhere. 
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In this letter we are going to consider a class of models where gravity interacts with a 
complex scalar field whose dynamics are dictated by a supergravity-Iike structure. More- 
over, we shall restrict ourselves to solving a set of first-order differential equations, the BPS 
equations, which can be deduced from the full, second-order, set of Einstein equations plus 
the equation of motion for the scalar field. Our goal is, given a Kahler potential and a 
superpotential for our scalar field, to find full numerical solutions for the field configuration 
and warp factors of the metric which localize gravity and give rise to an acceptable hier- 
archy between the Planck and electroweak scales in D=4. In addition, and based on our 
numerical results, we aim to establish general conditions for the existence of such trapping 
solutions. 

In the next section we present the models we want to study and the set of differential 
equations we want to solve. We also draw general conclusions about the structure of the 
solutions we are after. In section 3 we present models where the scalar field has canonical 
kinetic terms, solve the equations of motion and explain why the solutions obtained are not 
satisfactory. In section 4 we turn to models with a more complicated Kahler potential for 
which there are acceptable analytic solutions, and we explain their general characteristics. 
Finally, in section 5, we draw general results on the structure of these solutions and present 
the main conclusions. 



2. The model: gravity and a scalar field 



From now on we shall essentially follow the notation of ref. from where we also borrow 
their derivation of first-order BPS equations out of the second order Einstein plus field 
equations. We start off with an action in D-dimensions given by 



S = I d'^x VlffI ( - Kf.g-\da<f){db^) - V{<i^ 



ab I 



(2.1) 



where we work in units where 8ttGd = 1, fpi complex scalar fields, K^j = d'^ K / dcpidcpj 
is a Kahler metric, with K the Kahler potential. We are going to consider from now on 
one scalar field whose potential is given by the following expression 



V 



D 



{D - 2)K'^^W.^W.z - 2{D - 1)WW 



(2.2) 



where W{(j)) is a holomorphic function and plays the role of a superpotential. Moreover 

dW dK 



+ 



-W . 



(2.3) 



This form of the potential given by eq. (2.2) guarantees, at least for canonical Kahler 
metrics, vacuum stability with respect to small fluctuations of the AdSo geometry, see 



ref. H]. 

Once we have defined the matter content of our model, we can introduce the gravity 
sector. Since we are looking for string-like topological defects that could localize gravity, 
leading to our known four-dimensional world near the core of the string, we will fix from 
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now on D=6. A general metric consistent with four-dimensional Poincare invariance is 
given by 

ds^ = e^^^^'^'^r/^^dx'^dx^ + ie^^^^'^") {dzdz + dzdz) , (2.4) 

where = 1, ...4, z = + ix^, the signature of the D=4 metric is mostly plus and the 
two warp factors, A and B, are functions of the extra coordinates z, z. As we mentioned 
in the introduction, it is our intention to solve a reduced first-order version of the Einstein 



and scalar field equations, as derived in ref. |24], which is given by 



= e''''K^^W.^ , (2.5) 

d,A = -^e^*T^ , (2.6) 

d,N = K^d,^ + d,A = ie^* {2K'''^K^W.^ - W) , (2.7) 
where the function N is defined as 

N = B + ^ + iJ , (2.8) 



with J a real function. Note that equation (|2.6|) can be integrated once the system {\2,.5[ 
|2.7| ) has been solved. 

We are interested in axially symmetric solutions, therefore we will assume that, under a 
rotation z e*^z, the scalar field transforms as — > e*''^(/>, the Kahler potential is invariant 
and J gets shifted J ^ J + n9 for some integers q,n. Notice that this implies that W has 
to be a monomial, W = (f)^, and that 

n + pq = l. (2.9) 

In order to solve the system of equations give above we are going to parametrize the relevant 
functions as follows 

(j){z,z) = f{r)e''^^+''^^ , (2.10) 

^2A(z,,) ^ ^ (2.11) 

^N(z,z) ^ h(^ryj{r^)+ine+ieN ^ (2.12) 

where = zz, and we have implicitly assumed that J{z,z) = j{zz) + n9. For the time 
being we are going to set j{zz) = 0, although we shall return to this point later on. We 
want the functions /, g and h to be real and, for phenomenological purposes, we would like 



g to be positive-definite. Therefore, from eq. (2.6) we obtain a further constraint, namely 



eN+pO^ = 0, (2.13) 

modulo TT. We can always rotate away and, in such basis, Oj^i = 0. This is the choice we 
shall make from now on. 
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3. Canonical kinetic terms 



We are going to classify the models we want to study according to their Kahler potential. 
The simplest case is, then, that of considering canonical kinetic terms, i.e. K((j), (p) = (jxp. 
The scalar potential, eq. (|2.2|), is given by 



V{cl>,4>) = e*^2(</.0f-i[2(p+ - 5#] . (3.1) 

The system of BPS-equations becomes 

f + ^f = 2hr-\p + f), (3.2) 
ar r 

f + -h = 2h'mp + f)-l), (3.3) 
ar r 

% = -^9hr . (3.4) 

As we have already mentioned, the equation for g can be immediately integrated once 
the two first ones have been solved. Note that these equations are invariant under the 
transformations / — /, h (— l)'*'' h. Since the transverse metric is proportional h?{r), 
this function could be either positive or negative but, in any case, has to be different from 
zero. 

The condition for trapping gravity is given by 

dx^dx^g^^^M ^ f dr e^^+^B^ = ( dr g{r)h^ {r)e- < oo . (3.5) 



If we assume that the solution is regular everywhere (including the origin) the condition 
for a finite Planck mass depends on the behaviour of the configuration for large r values. 
It is easy to check that a possible consistent ansatz in the large r region is 

/(r)|™ ~ /o , (3.6) 

h{r)\r-,oo ~ — , (3.7) 
r 

with /o, ho constants. From eq. (|3.2| ) we can calculate 

r2-p 

We can also calculate g using this ansatz. The result is 

5(r)|™~<7or-2'^°^o . (3.9) 



^0 = • (3.8) 



Then, the integral given by eq. ( |3.5| ) is convergent when 

/io/o^>0 , (3.10) 

i.e. when 



^ >0 (3.11) 
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is fulfilled. 

As mentioned above, the validity of this ansatz relies on the existence of a real solution 
to the algebraic equation for /q. Combining eqs. (|3.2| , ^ ) we get 



1 



(3.12) 



In the preceding discussion, the value of p, the degree of the superpotential term, was 
an arbitrary integer. From eq. ( 3.121) we deduce that a real value for /o implies p < 0. 
However, a non singular scalar potential implies p > 1 or p = 0; therefore we will fixp = 0. 
Note that the trapping condition, eq. ( |3.11| ), for that case translates into g > 0. 

After studying the necessary conditions to have gravity trapping in the large r region, 
we should check the hypothesis of regularity at r = 0. Once we fix = 1 it turns out that 
there are no regular, consistent, boundary conditions at the origin which are compatible 
with this choice. This can be also inferred from the shape of the potential 



5) 



(3.13) 



shown in fig. |T[ There we can see that, at least up to the minimum, the value of the 
potential energy is negative. This is a crucial fact that prevents the model from working 
and which we shall explain in detail in section 5. 




Figure 1: Scalar potential, as given in eq. (3.13), as a function of 



4. Non-canonical kinetic terms 



Motivated by the previous discussion, we have generalized the analysis to the case where 
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the Kahler potential is not canonical. The BPS equations for this case are 

f + ^/ = W'-'i^, (4.1) 



dr r k' + f^k" ' 

dr r k' + f^k 



+ -h = h'Fi2k'J;^^ - 1) , (4.2) 



I = , (4.3) 

where k{4>(j)) = K{(j), 0) and ' denotes a derivative with respect to (fxp. The corresponding 
equations (|3.2| - |3.4[) for the canonical case are recovered by taking k' = l,k" = 0. 

The equivalent to equations ( |3.8| ) and ( f3.12 ) for the asymptotic ansatz are given by 



where A;^ = fc'(/2),A:^' = A;"(/o2). 

For a regular Kahler potential, the trapping condition is the same as the one for the 
canonical case, eq. ( |3.1C1| ), namely /iq/q > which, in terms of /o, reads 

Once we have determined all possible phenomenological constraints on the model, we should 
try to look for solutions of the BPS equations that fulfill them. Although we have not 
performed a systematic analysis, we have tried several simple analytical ansatze for the 
Kahler potential as a function of which, together with our choice for W, define the 
particle content. Within this context it is not difficult to find examples where all constraints 
are satisfied. As we will see in the next section, once we fix (p, q), the boundary conditions 
at the origin are essentially unique (modulo a rescaling of the coordinate z). Then we just 
have to evolve the different functions from these regular boundary conditions at r = to 
the large r regime. The problem is that, in most cases, the solutions do not have a good 
asymptotic behaviour, i.e. the condition for localizing gravity in 4D is no longer fulfilled. 

4.1 A particular model 

A simple working example is described by 

W(4>) = * , (4.8) 
which corresponds to choosing p = 1, and gives the following scalar potential 
V(, I) - .Jt^) ^ - ^^^^ + 11^^'^''^' - 7100'^^'^=^ + 20000</>^<^^ 

^^'^^ 1-1000202 • ^ • ^ 
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It is worth pointing out that the metric derived from this choice of Kahler potential becomes 
singular at = 0.1 and is well defined only for smaller values of The potential is 
positive for small \(j)\ values and then decreases, entering a region where V is negative. It 
has a minimum at \(f)\'^ ~ 0.097, close to the barrier due to the above-mentioned singularity 
in the Kahler potential^. All these features are shown in fig. |2|. Let us now try to solve 

25 
20 
15 

> 5 
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Figure 2: Plot of the scalar potential as a function of \(f>\^ for the model described in the text. 

the BPS equations. Our ansatz for the scalar field and one of the warp factors is 

<t> = f{ry' , (4.10) 
e^ = h{r) . (4.11) 

Note that this choice corresponds to g = 1 and n = (i.e. J does not depend on 9). By 
solving the equations for /o and /iq, we get 

/o-f,/..-^. (4.12) 

As we have mentioned before, /q = 0.08 < 0.1 and, therefore, the asymptotic value for |(/)| 
falls in the region where the model is well defined. Notice that V < Q for |(/)| = /o, which 
means that the potential plays the role of a negative cosmological constant in the external 
region. 



^This singularity does not pose a problem since, as we will see, the string does not spread beyond 
such value. In fact we have found other, more involved, models where the Kahler metric is well defined 
everywhere. For example, one can take the Kahler potential presented here, perform a series expansion 
around the value \<f)\'^ = 0, and keep as many terms as necessary in order to reproduce the scalar potential 
shape given in fig. |^ in the region of interest. Such potentials are perfectly well behaved, though one must 
keep the first eight terms of the Kahler expansion in order to reproduce all the features of the model. 
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Concerning the string core, we have found that there is only one consistent set of 
regular boundary conditions^. This is given by / ~ r, /i ~ const. Since it is possible to 
rescale z by redefining A, B we can fix, at lowest order, f = r. We then get 



f{r) =r + 0{r') , (4.13) 

2 

/i(r) = 1 + y + 0(r3) , (4.14) 

Mr) = -^+0{r^) , (4.15) 

B{r) = 0{r^) . (4.16) 
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Figure 3: Plot of f{r)/fo as a function of r. 

Once we have fixed these boundary conditions, we just integrate the differential equa- 
tions. The result is shown in figs. ^ ^. As we anticipated, both functions approach their 
asymptotic values in the large r region. To illustrate this fact, we have plotted f{r)/fo in 
fig. ^ and, in fig. ^, we have added the theoretical asymptotic prediction for h{r), namely 
ho/r. Notice that the string core deduced from the scalar field stabilization, /(r), is 0(10) 
smaller than the one given by the metric variation, h{r). 

Using eqs. ( |2.8D , ( |2.1l| ) and ( ^.12| ) we can evaluate the metric factors e^"^'^^^ e^^*^''^. 
These are shown in fig. |5[ For large r values 

e2^(^) ~ 2r-'^o^o = , (4.17) 

^ These boundary conditions apply only to the case q = 1, however we have found that there exists only 
one solution for any q > 0, which is nothing but the one presented here expressed in terms of the variable 
z = u". 
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Figure 4: Plot of h{r) and ho/r as a function of r. 




50 100 150 200 



r 

Figure 5: Plot of the two warp factors, e'^'''-' and e'^'''-' as a function of 



Let us now compute the energy of the configuration described above. We will use 
the notation and results given in [^]. For static configurations, the energy stored in the 
transverse plane is defined as 

E = Egra.v + -Ekin + -Epot (4.19) 

and is equal to the negative of the action, as given by eq. (|2.1D. In ref. [24| is shown that 



it is possible to write the energy density as a sum of quadratic terms, each vanishing when 
the BPS equations are fulfilled, plus some total derivatives Si,S2,S3 given by 

51 = 5[e^^(8 dA + 2dN- dK)] + {d^B ,N ^N*), (4.20) 

52 = d[e^^{K^dcl) - K^d4>)] + (d^B ,(j)^4>) , (4.21) 

53 = 45[e^^+^ W] + 4a[e^^+^* W] . (4.22) 

Then the energy of a BPS configuration is given by 

E = j dx^cix^(Si + S2 + S3). (4.23) 

We have found that it is still possible to simplify this expression, using the BPS equations, 
once again. We get 

E = J dx^ dx^ {2d[e^'^{BA + B{N - N*))] + {d^B ,N ^ N*)) . (4.24) 
Using Stokes' theorem 

/ dz f\dz{dzV2 — d^Vz) = f {dzvz + dzvz), (4.25) 

J R Jc=SFt^ 

with 

Vz = i 2e'^^dz{A + N - N*) , (4.26) 
we can calculate the energy stored in a region enclosed by a contour C 

Ec = -i j (2e^'^dz{A + N* - N)dz - 2e^^d-z{A + N - N*)dz) . (4.27) 

Let us now aply this formula to an axially symmetric solution. The first thing to notice is 
that the contribution from N cancels when this function is real, as happens to be the case 
in our example^. In that situation the energy can be seen as uniquely dependent on the 
warp factor. 
Using that 

zdze^^^^'^'^ =rg{r)g'{r), (4.28) 



^In general this term will vanish if the two real functions {h,j) that parametrize A'^ have no singularities, 
as expected in regular solutions. 
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we can calculate the energy E(r) stored in a circle of radius r around the center of the 
string 

E{r) =27rrg{r)g'{r) . (4.29) 

We have verified this equation numerically. Taking the limit r ^ co we find that Etotai = 0. 

We should stress that this energy is zero for any solution of that traps gravity. This 
follows from the regularity of the metric at the origin, r = 0, and from the inverse power 
decay of the warp factor, A{r). 




2 4 6 8 10 12 14 16 



r 



Figure 6: Plot of the energy enclosed in a circle of radius r as a function of r. The total energy 
E{r) is given by the thick black line, whereas its three different terms (i?grav(f) -blue-, Ekinir) 
-red- and £'pot('") -green-) are also shown. 

We illustrate the previous results in fig. where we have plotted E{r), calculated 
by numerical integration, as a function of r (thick black line). We have also plotted the 
three different contributions to the energy, coming from the integration of the curvature 
(blue solid line), kinetic (red) and potential (green) energy densities as given in ( 4.2C| ). 



A remarkable fact is that the total (integrated) potential energy goes to zero. This can 
be easily understood by calculating the energy of a family of configurations that are just 
scaling deformations of the solution, and using the fact that our defect spreads along two 
spatial dimensions. The vanishing of the total potential energy sets stringent bounds on 
the class of models of gravity plus scalar field that could work. In particular, models with 
canonical kinetic terms, as presented in section 3, are ruled out, as it is now obvious from 
fig. [l|, where the potential happens to be negative in the region where the string would 
spread. 



5. Discussion 

In this section we would like to elaborate further on the main results achieved. In order 
to do that, let us present first the geometry that results of our solution. To see it more 
straightforwardly, it is useful to realize that there are other parametrizations of metrics 
in six dimensions with axial symmetry. A popular one is given by a normalized radial 
contribution to the metric, i.e 

ds^ = a^{p) r]^^dx>'dx'' + dp'^ + j'^{p) dO^ . (5.1) 

We have expressed our solutions in terms of these coordinates. The change of variables is 
regular everywhere. The shape of the cr{p) and 7(p) factors is drawn in fig. |7[ Far from 

1 
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P 

Figure 7: a{p) and j{p)/Rq as a function of p for the solution presented in the previous section, 
the core, the a warp factor decreases exponentially 

a\p)r^e~P\ (5.2) 

with 

and 7 goes to a constant given by 

i?o = |/io|e-'=[^ol/2 

The space factorizes as AdS^ x Si and, therefore, the curvature in this region is constant. 
This geometry, with a cylindrical transverse space and an exponential warp factor, was 




(5.3) 
(5.4) 
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proposed by Gregory Q as the one allowing for nonsingular global string compactifica- 
tions in D=6, and also by Olasagasti and Vilenkin [^] in their generalization to higher 
dimensional global defects. 

Inside the core, 7(p) ~ p and, consequently, there is no deficit angle. This follows from 
our choice of boundary conditions at the origin, as can be seen using that p ~ r for small 
values of p. 

Our solution is regular everywhere and interpolates smoothly between the origin and 
the asymptotic region, providing a cigar-like, complete geometry. As we mentioned before, 
this kind of geometry was proposed by Gregory in ref. Q. She considered a global string 
given by a potential with Vmm = and a negative cosmological constant and argued that, 
for some values of the parameters, it should be possible to match the two regions building 
a nonsingular metric. We have presented here an example based on this idea, although 
there is one difference which is important enough to be pointed out. In our example there 
is no cosmological constant. The required negative contribution to the energy density far 
from the core is provided by the potential and, therefore, depends on the asymptotic value 
for \(f)\. This value is not necessarily, as one could naively think, the one that minimizes 
the potential. It is dynamically fixed by the field equations which implies that, in general, 
we can not neglect the feedback of the metric on the stabilization of the modulus of the 
scalar field. In fig. we plot the gravitational, kinetic and potential energy densities, as 
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Figure 8: Plot of the different energy densities as a function of p. As in fig. ^ here the total 
energy density is given by the thick black line, whereas its different contributions are given by the 
blue (curvature), red (kinetic) and green (potential) lines. 

well as their sum. The pattern shown in the figure is quite general. As explained before, 
far from the core the value of the curvature is fixed by the scalar -winding plus potential- 
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energy, A,^, through the AdS relation 

i? = ^A. (5.5) 

whith d = 5. We can check that the negative potential energy density dominates over the 
kinetic one in this asymptotic region. Since the energy coming from the integration of the 
curvature term has to cancel the total kinetic energy, the curvature must be positive in 
the inner region, as shown in fig. |^. The value of p at which the scalar curvature vanishes 
{p = 2) is of the same order of magnitude than the one where V also does (p = 1.3). 

The previous discussion leads to the next issue we want to address, namely the phe- 
nomenological applications of these models. From the features of the energy configurations 
cited above, it turns out that we can observe how the relevant parameters of the model 
(i.e. the two warp factors and the modulus of the scalar field) adjust themselves to their 
asymptotic values, a, Rq and /o, in order to fulfill them. In fact, we can relate the asymp- 
totic value of a to the value of the scalar potential at /o, ^[/q], as we have written in 
eq. (H. 

This means, among other things, that obtaining the right hierarchy between the D=6 
and D=4 Planck masses (mainly parameterized by a), will crucially depend on the dy- 
namics of the scalar field, namely its Kahler potential and its superpotential. Indeed the 
trapping condition eq. ( 3.1C| ) has to be verified in order to have a finite four-dimensional 



reduced Planck mass. Also the transverse volume has to be large enough to generate the 
appropriate value for the Planck mass from the six-dimensional parameters (these are the 
higher-dimensional scale of gravity, M, and the parameters upon which W and K depend). 
In the case we have presented this volume (see eq. ( |3.5| )) is O(IOO), therefore a gravity scale 
hierarchy between D=6 and D=4 is not generated. 

We have not yet estimated the correction to Newton's law in this kind of scenarios. 
In general we expect a mixture of effects that typically appear in the case of one extra 
transverse dimension plus some Kaluza-Klein ones modulated by the presence of the warp 
factor. Since we have a complete description of the metric, one could perform, at least 
numerically, a complete analysis within a realistic model. Meanwhile, we can use the 
results of two scenarios that share some of the features of our model. In ref. [^], a single, 
positive tension four-brane embedded in a six-dimensional AdS space is considered. One 
transverse dimension is compactified and the other is taken to be an infinite one. They 
consider axially symmetric, constant energy momentum tensors. These bulk and brane 
tensors are allowed to be inhomogeneous and include, in particular, the typical pattern of 
asymptotic string-like defects [10|. The flatness of the brane and the specific cylindrical 



geometry are obtained by several fine tunings among the bulk and brane energy momentum 
tensor components. The correction to Newton's law in this scenario is shown to be small. 
Another particular case, providing this kind of energy momentum tensor, is obtained by 
considering an Abelian magnetic field [^, The radius of the cylinder is fixed by the 
magnetic charge and by the bulk cosmological constant. 

We want to stress that in these models there is a core where the metric is unknown 
since the asymptotic solution is not valid there. On the other hand, in order to evaluate 
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the corrections to Newton's law, it is either necessary to make some assumptions on the 
form of the metric near p = 0, or one has to work in the thin-core hmit. Since this latter 
approach is not always justified, our solution could provide a model to explore general thick 
core corrections to these approximations. 

There are other issues that one has to study when the core of the defect is replaced 
by a brane with some tension. In D=5, for example, one is almost forced to include 
negative tension branes to generate decreasing warp factors. This is dangerous, since the 
energy conditions will, in general, be violated. Also one has impose energy momentum 
conservation by hand, getting relations between brane tensions and metric discontinuities. 
Since we are dealing with a real topological defect, our model does not suffer this kind of 
anomalies. The metric is regular everywhere and the energy momentum tensor is conserved 
- since it is derived from a Lagrangian density -. Also the null energy condition 

Tapv'^v'^ > , (5.6) 

with a null vector, is automatically fulfilled since 

To^pv^vf" = 2|^^|2 K^-^ {id<j>)id4>) + (5.7) 

is always positive. 

In conclusion we have presented here a complete numerical solution to the system 
of Einstein and field equations of a six-dimensional setup with gravity and a scalar field. 
This was achieved by solving a particular set of first order differential equations, the BPS 
equations, whose solutions also satisfy the more general second order system. Using analytic 
ansatze on the behaviour of the field and warp factors of the metric around the origin 
of the defect and in the aymptotic region, we have been able to classify the solutions 
obtained and determine what sort of dynamics (i.e. Kahler potential and superpotential) 
must the scalar field obey in order to have gravity trapping around the four-dimensional 
defect. The geometry of the phenomenologically viable solutions is cigar-shaped and those 
configurations have a total vanishing energy. 
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